Protein Dynamics Research Articles

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics.

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Conformational Dynamics in Corynebacterium glutamicum Diaminopimelate Epimerase: Insights from Ligand Parameterization, Atomistic Simulation, and Markov State Modeling.

The microbial enzyme diaminopimelate epimerase (DapF), a vital enzyme in the lysine biosynthetic pathway, catalyzes the conversion of L, L-diaminopimelate (L, L-DAP) to D, L-diaminopimelate (D, L-DAP) using a catalytic cysteine dyad with one cysteine in thiol state and another in thiolate. Under oxidizing conditions, the catalytic cysteines of apo DapF form a disulfide bond that alters the structure and function of DapF. Given its potential as a target for antimicrobial resistance treatments, understanding DapF's functional dynamics is imperative. In the present work, we employ microsecond-scale all-atom molecular dynamics simulations of product-bound DapF and apo-DapF under oxidized and reduced conditions. We employ a polarized charge model for the ligand and the active site residues, which was necessary to preserve the electrostatic environment in the active site and retain the ligand in the active site. The product-bound DapF and apo-DapF in oxidized and reduced conditions exhibit a closed, semi-open, and open conformation, respectively, as identified using the internal coordinates of the dimeric enzyme and the principal component analysis. The conformational switch is guided by the dynamic catalytic (DC) loop, loop II, and loop III movements in the active site. The time scale of the close-to-open conformational transition is estimated to be 0.8 μs through Markov state modeling (MSM) and transition path theory (TPT). The present study explains the role of various active site residues and loops in ligand binding and protein dynamics in the DapF enzyme under different redox conditions. Such information will be helpful in future inhibitor design studies targeting the DapF enzyme.

The Applications of Molecular Dynamics Simulation in Studying Protein Structure and Dynamics

The Applications of Molecular Dynamics Simulation in Studying Protein Structure and Dynamics

Dynamics of the left ventricular ejection fraction after revascularization in patients with heart failure with preserved ejection fraction, association with the TRAIL protein

Aim. To study the association of dynamics of the ejection fraction (EF) of the left ventricle (LV) with the development of adverse cardiovascular events within 12 months after revascularization in patients with chronic heart failure with preserved LV EF (HFpEF), and to determine the value of TNF-related apoptosis-inducing ligand (TRAIL) in predicting changes of LVEF. Materials and methods. 52 patients with HFpEF hospitalized for coronary artery bypass grafting (CABG) were included in the prospective study. The levels of plasma NTproBNP (before CABG) and TRAIL protein (before CABG and 10 days after CABG) were determined. Echocardiography was performed for these patients at the time of enrollment in the study and after 12 months of follow-up. The patients were divided into 2 groups: group 1 (n=23) included patients with increased LVEF, group 2 (n=29) – with unchanged or reduced LVEF. The development of a combined endpoint (cardiovascular death, decompensated of HF, repeat unplanned revascularization) was studied. Results. LVEF before CABG was comparable in the studied groups. The frequency of the combined endpoint was 0 in the first group and 17.2% in the second group (p=0.02). There was an inverse correlation between change of LVEF one year after revascularization and the dynamics of the TRAIL protein in the perioperative period (r=-0.45, p=0.049). The area under the ROC for perioperative TRAIL changes and changes in LVEF during a follow-up was 0.79 (95% confidence interval 0.6–0.985; p=0.018). Conclusion. The increase of LVEF 12 months after CABG in patients with HFpEF is associated with a rarer development of cardiovascular events. The changes of LVEF after CABG can be predicted based on the analysis of the perioperative changes of the TRAIL protein.

Structural Characterization of Monoclonal Antibodies and Epitope Mapping by FFAP Footprinting.

Covalent labeling in combination with mass spectrometry is a powerful approach used in structural biology to study protein structures, interactions, and dynamics. Recently, the toolbox of covalent labeling techniques has been expanded with fast fluoroalkylation of proteins (FFAP). FFAP is a novel radical labeling method that utilizes fluoroalkyl radicals generated from hypervalent Togni reagents for targeting aromatic residues. This report further demonstrates the benefits of FFAP as a new method for structural characterization of therapeutic antibodies and interaction interfaces of antigen-antibody complexes. The results obtained from human trastuzumab and its complex with human epidermal growth factor receptor 2 (HER2) correlate well with previously published structural data and demonstrate the potential of FFAP in structural biology.

Dynamic changes in proteomics during the response of Arabidopsis seedlings to heat stress.

Abstract The threat of global warming to plant survival, along with its adverse effects on growth and agricultural productivity, underscores the urgency of developing thermotolerant crops. Achieving this goal necessitates a thorough comprehension of plant responses to heat stress at the molecular level. To address this, we conducted an investigation into proteome dynamics in Arabidopsis thaliana seedlings exposed to moderate heat stress (30°C). Employing a novel approach combining 15N-stable isotope labeling and the ProteinTurnover algorithm facilitated a detailed examination of proteomic changes across various cellular fractions. Our study unveiled significant alterations in the turnover rates of 571 proteins, with a median increase of 1.4-fold, indicative of heightened protein dynamics in response to heat stress. Notably, soluble proteins in the roots exhibited comparatively minor changes, suggesting tissue-specific adaptive mechanisms. The analysis also revealed substantial turnover variations in proteins associated with redox signaling, stress response, and metabolism, underscoring the intricate nature of the response network. Conversely, proteins involved in carbohydrate metabolism and mitochondrial ATP synthesis demonstrated minimal turnover fluctuations, highlighting their inherent stability. This comprehensive assessment sheds light on the proteomic adaptations of Arabidopsis seedlings to moderate heat stress, elucidating the delicate balance between proteome stability and adaptability. These insights contribute to our understanding of plant thermal resilience and offer valuable support for the development of crops endowed with enhanced thermotolerance.

Interpretation of Hydrogen/Deuterium Exchange Mass Spectrometry.

This paper sheds light on the meaning of hydrogen/deuterium exchange-mass spectrometry (HDX-MS) data. HDX-MS data provide not structural information but dynamic information on an analyte protein. First, the reaction mechanism of backbone amide HDX reaction is considered and the correlation between the parameters from an X-ray crystal structure and the protection factors of HDX reactions of cytochrome c is evaluated. The presence of H-bonds in a protein structure has a strong influence on HDX rates which represent protein dynamics, while the solvent accessibility only weakly affects the HDX rates. Second, the energy diagrams of the HDX reaction at each residue in the presence and absence of perturbation are described. Whereas the free energy change upon mutation can be directly measured by the HDX rates, the free energy change upon ligand binding may be complicated due to the presence of unbound analyte protein in the protein-ligand mixture. Third, the meanings of HDX and other biophysical techniques are explained using a hypothetical protein folding well. The shape of the protein folding well describes the protein dynamics and provides Boltzmann distribution of open and closed states which yield HDX protection factors, while a protein's crystal structure represents a snapshot near the bottom of the well. All biophysical data should be consistent yet provide different information because they monitor different parts of the same protein folding well.

Impact of N-Glycosylation on Protein Structure and Dynamics Linked to Enzymatic C-H Activation in the M. oryzae Lipoxygenase.

Lipoxygenases (LOXs) from pathogenic fungi are potential therapeutic targets for defense against plant and select human diseases. In contrast to the canonical LOXs in plants and animals, fungal LOXs are unique in having appended N-linked glycans. Such important post-translational modifications (PTMs) endow proteins with altered structure, stability, and/or function. In this study, we present the structural and functional outcomes of removing or altering these surface carbohydrates on the LOX from the devastating rice blast fungus, M. oryzae, MoLOX. Alteration of the PTMs did notinfluence the active site enzyme-substrate ground state structures as visualized by electron-nuclear double resonance (ENDOR) spectroscopy. However, removal of the eight N-linked glycans by asparagine-to-glutamine mutagenesis nonetheless led to a change in substrate selectivity and an elevated activation energy for the reaction with substrate linoleic acid, as determined by kinetic measurements. Comparative hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis of wild-type and Asn-to-Gln MoLOX variants revealed a regionally defined impact on the dynamics of the arched helix that covers the active site. Guided by these HDX results, a single glycan sequon knockout was generated at position 72, and its comparative substrate selectivity from kinetics nearly matched that of the Asn-to-Gln variant. The cumulative data from model glyco-enzyme MoLOX showcase how the presence, alteration, or removal of even a single N-linked glycan can influence the structural integrity and dynamics of the protein that are linked to an enzyme's catalytic proficiency, while indicating that extensive glycosylation protects the enzyme during pathogenesis by protecting it from protease degradation.

Secondary Sites of the C-type Lectin-Like Fold.

C-type lectins are a large superfamily of proteins involved in a multitude of biological processes. In particular, their involvement in immunity and homeostasis has rendered them attractive targets for diverse therapeutic interventions. They share a characteristic C-type lectin-like domain whose adaptability enables them to bind a broad spectrum of ligands beyond the originally defined canonical Ca2+-dependent carbohydrate binding. Together with variable domain architecture and high-level conformational plasticity, this enables C-type lectins to meet diverse functional demands. Secondary sites provide another layer of regulation and are often intricately linked to functional diversity. Located remote from the canonical primary binding site, secondary sites can accommodate ligands with other physicochemical properties and alter protein dynamics, thus enhancing selectivity and enabling fine-tuning of the biological response. In this review, we outline the structural determinants allowing C-type lectins to perform a large variety of tasks and to accommodate the ligands associated with it. Using the six well-characterized Ca2+-dependent and Ca2+-independent C-type lectin receptors DC-SIGN, langerin, MGL, dectin-1, CLEC-2 and NKG2D as examples, we focus on the characteristics of non-canonical interactions and secondary sites and their potential use in drug discovery endeavors.

Brain chains as topological signatures for Alzheimer’s disease

AbstractTopology is providing new insights for neuroscience. For instance, graphs, simplicial complexes, directed graphs, flag complexes, persistent homology and convex covers have been used to study functional brain networks, synaptic connectivity, and hippocampal place cell codes. We propose a topological framework to study the evolution of Alzheimer’s disease, the most common neurodegenerative disease. The modeling of this disease starts with the representation of the brain connectivity as a graph and the seeding of a toxic protein in a specific region represented by a vertex. Over time, the accumulation of toxic proteins at vertices and their propagation along edges are modeled by a dynamical system on this graph. These dynamics provide an order on the edges of the graph according to the damage created by high concentrations of proteins. This sequence of edges defines a filtration of the graph. We consider different filtrations given by different disease seeding locations. To study these filtrations we propose a new combinatorial and topological method. A filtration defines a maximal chain in the partially ordered set of spanning subgraphs ordered by inclusion. To identify similar graphs, and define a topological signature, we quotient this poset by graph homotopy equivalence, which gives maximal chains in a smaller poset. We provide an algorithm to compute this direct quotient without computing all subgraphs and then propose bounds on the total number of graphs up to homotopy equivalence. To compare the maximal chains generated by this method, we extend Kendall’s $$d_K$$ d K metric for permutations to more general graded posets and establish bounds for this metric. We then demonstrate the utility of this framework on actual brain graphs by studying the dynamics of tau proteins on the structural connectome. We show that the proposed topological brain chain equivalence classes distinguish different simulated subtypes of Alzheimer’s disease.

Interplay Between Substrate Polarity and Protein Dynamics in Evolved Kemp Eliminases

Extensive molecular dynamics (MD) simulations on designed and evolved enzymes for the Kemp Elimination were performed. Over thirty different systems including combinations of three protein scaffolds used in computational design (KE70, KE59 and HG3), several directed evolution variants and ligands (5‐nitro‐benzoisoxazole, and 5‐nitro‐benzotriazole), were screened computationally. The study was focused mostly on computational and evolved variants for which X‐ray structures are available. MD simulations were used to monitor the organization of the active sites over time periods of 100 ns. A general trend was observed in the simulations: the optimal organization of active site residues, corresponding to the most active catalysts, is only maintained with polarized, TS‐like substrates, reinforcing the crucial role of protein dynamics and electrostatic preorganization for efficient biocatalysis. Implications for enzyme design protocols are discussed.

Dynamics of Single-Cell Protein Covariation during Epithelial-Mesenchymal Transition.

Physiological processes, such as the epithelial-mesenchymal transition (EMT), are mediated by changes in protein interactions. These changes may be better reflected in protein covariation within a cellular cluster than in the temporal dynamics of cluster-average protein abundance. To explore this possibility, we quantified proteins in single human cells undergoing EMT. Covariation analysis of the data revealed that functionally coherent protein clusters dynamically changed their protein-protein correlations without concomitant changes in the cluster-average protein abundance. These dynamics of protein-protein correlations were monotonic in time and delineated protein modules functioning in actin cytoskeleton organization, energy metabolism, and protein transport. These protein modules are defined by protein covariation within the same time point and cluster and, thus, reflect biological regulation masked by the cluster-average protein dynamics. Thus, protein correlation dynamics across single cells offers a window into protein regulation during physiological transitions.

Structure and Mechanisms of Assembly-Line Polyketide Synthases.

Three decades of studies on the multifunctional 6-deoxyerythronolide B synthase have laid a foundation for understanding the chemistry and evolution of polyketide antibiotic biosynthesis by a large family of versatile enzymatic assembly lines. Recent progress in applying chemical and structural biology tools to this prototypical assembly-line polyketide synthase (PKS) and related systems has highlighted several features of their catalytic cycles and associated protein dynamics. There is compelling evidence that multiple mechanisms have evolved in this enzyme family to channel growing polyketide chains along uniquely defined sequences of 10-100 active sites, each of which is used only once in the overall catalytic cycle of an assembly-line PKS. Looking forward, one anticipates major advances in our understanding of the mechanisms by which the free energy of a repetitive Claisen-like reaction is harnessed to guide the growing polyketide chain along the assembly line in a manner that is kinetically robust yet evolutionarily adaptable.

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation.

Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.

Maximum entropy determination of mammalian proteome dynamics

Full understanding of proteostasis and energy utilization in cells will require knowledge of the fraction of cell proteins being degraded with different half-lives and their rates of synthesis. We therefore developed a method to determine such information that combines mathematical analysis of protein degradation kinetics obtained in pulse-chase experiments with Bayesian data fitting using the maximum entropy principle. This approach will enable rapid analyses of whole-cell protein dynamics in different cell types, physiological states, and neurodegenerative disease. Using it, we obtained surprising insights about protein stabilities in cultured cells normally and upon activation of proteolysis by mTOR inhibition and increasing cAMP or cGMP. It revealed that >90% of protein content in dividing mammalian cell lines is long-lived, with half-lives of 24 to 200 h, and therefore comprises much of the proteins in daughter cells. The well-studied short-lived proteins (half-lives < 10 h) together comprise <2% of cell protein mass, but surprisingly account for 10 to 20% of measurable newly synthesized protein mass. Evolution thus appears to have minimized intracellular proteolysis except to rapidly eliminate misfolded and regulatory proteins.

Fatty acid-binding protein 5 is a functional biomarker and indicator of ferroptosis in cerebral hypoxia

The progression of human degenerative and hypoxic/ischemic diseases is accompanied by widespread cell death. One death process linking iron-catalyzed reactive species with lipid peroxidation is ferroptosis, which shows hallmarks of both programmed and necrotic death in vitro. While evidence of ferroptosis in neurodegenerative disease is indicated by iron accumulation and involvement of lipids, a stable marker for ferroptosis has not been identified. Its prevalence is thus undetermined in human pathophysiology, impeding recognition of disease areas and clinical investigations with candidate drugs. Here, we identified ferroptosis marker antigens by analyzing surface protein dynamics and discovered a single protein, Fatty Acid-Binding Protein 5 (FABP5), which was stabilized at the cell surface and specifically elevated in ferroptotic cell death. Ectopic expression and lipidomics assays demonstrated that FABP5 drives redistribution of redox-sensitive lipids and ferroptosis sensitivity in a positive-feedback loop, indicating a role as a functional biomarker. Notably, immunodetection of FABP5 in mouse stroke penumbra and in hypoxic postmortem patients was distinctly associated with hypoxically damaged neurons. Retrospective cell death characterized here by the novel ferroptosis biomarker FABP5 thus provides first evidence for a long-hypothesized intrinsic ferroptosis in hypoxia and inaugurates a means for pathological detection of ferroptosis in tissue.

Pooled multicolour tagging for visualizing subcellular protein dynamics.

Imaging-based methods are widely used for studying the subcellular localization of proteins in living cells. While routine for individual proteins, global monitoring of protein dynamics following perturbation typically relies on arrayed panels of fluorescently tagged cell lines, limiting throughput and scalability. Here, we describe a strategy that combines high-throughput microscopy, computer vision and machine learning to detect perturbation-induced changes in multicolour tagged visual proteomics cell (vpCell) pools. We use genome-wide and cancer-focused intron-targeting sgRNA libraries to generate vpCell pools and a large, arrayed collection of clones each expressing two different endogenously tagged fluorescent proteins. Individual clones can be identified in vpCell pools by image analysis using the localization patterns and expression level of the tagged proteins as visual barcodes, enabling simultaneous live-cell monitoring of large sets of proteins. To demonstrate broad applicability and scale, we test the effects of antiproliferative compounds on a pool with cancer-related proteins, on which we identify widespread protein localization changes and new inhibitors of the nuclear import/export machinery. The time-resolved characterization of changes in subcellular localization and abundance of proteins upon perturbation in a pooled format highlights the power of the vpCell approach for drug discovery and mechanism-of-action studies.

Parkinsonism Sac domain mutation in Synaptojanin-1 affects ciliary properties in iPSC-derived dopaminergic neurons

Synaptojanin-1 (SJ1) is a major neuronal-enriched PI(4, 5)P2 4- and 5-phosphatase implicated in the shedding of endocytic factors during endocytosis. A mutation (R258Q) that impairs selectively its 4-phosphatase activity causes Parkinsonism in humans and neurological defects in mice (SJ1RQKI mice). Studies of these mice showed, besides an abnormal assembly state of endocytic factors at synapses, the presence of dystrophic nerve terminals selectively in a subset of nigro-striatal dopamine (DA)-ergic axons, suggesting a special lability of DA neurons to the impairment of SJ1 function. Here we have further investigated the impact of SJ1 on DA neurons using iPSC-derived SJ1 KO and SJ1RQKI DA neurons and their isogenic controls. In addition to the expected enhanced clustering of endocytic factors in nerve terminals, we observed in both SJ1 mutant neuronal lines increased cilia length. Further analysis of cilia of SJ1RQDA neurons revealed abnormal accumulation of the Ca2+ channel Cav1.3 and of ubiquitin chains, suggesting a defect in the clearing of ubiquitinated proteins at the ciliary base, where a focal concentration of SJ1 was observed. We suggest that SJ1 may contribute to the control of ciliary protein dynamics in DA neurons, with implications on cilia-mediated signaling.

Hydrophobic Clusters Regulate Surface Hydration Dynamics of Bacillus subtilis Lipase A.

The surface hydration diffusivity of Bacillus subtilis Lipase A (BSLA) has been characterized by low-field Overhauser dynamic nuclear polarization (ODNP) relaxometry using a series of spin-labeled constructs. Sites for spin-label incorporation were previously designed via an atomistic computational approach that screened for surface exposure, reflective of the surface hydration comparable to other proteins studied by this method, as well as minimal impact on protein function, dynamics, and structure of BSLA by excluding any surface site that participated in greater than 30% occupancy of a hydrogen bonding network within BSLA. Experimental ODNP relaxometry coupling factor results verify the overall surface hydration behavior for these BSLA spin-labeled sites similar to other globular proteins. Here, by plotting the ODNP parameters of relative diffusive water versus the relative bound water, we introduce an effective "phase-space" analysis, which provides a facile visual comparison of the ODNP parameters of various biomolecular systems studied to date. We find notable differences when comparing BSLA to other systems, as well as when comparing different clusters on the surface of BSLA. Specifically, we find a grouping of sites that correspond to the spin-label surface location within the two main hydrophobic core clusters of the branched aliphatic amino acids isoleucine, leucine, and valine cores observed in the BSLA crystal structure. The results imply that hydrophobic clustering may dictate local surface hydration properties, perhaps through modulation of protein conformations and samplings of the unfolded states, providing insights into how the dynamics of the hydration shell is coupled to protein motion and fluctuations.

USING FÖRSTER RESONANCE ENERGY TRANSFER (FRET) TO UNDERSTAND THE UBIQUITINATION LANDSCAPE.

Förster resonance energy transfer (FRET) is a fluorescence technique that allows quantitative measurement of protein interactions, kinetics and dynamics. This review covers the use of FRET to study the structures and mechanisms of ubiquitination and related proteins. We survey FRET assays that have been developed where donor and acceptor fluorophores are placed on E1, E2 or E3 enzymes and ubiquitin (Ub) to monitorsteady-stateand real-time transfer of Ub through the ubiquitination cascade. Specialized FRET probesplaced onUb andUb-likeproteins have been developed to monitor Ub removal by deubiquitinating enzymes (DUBs)that result in a loss of aFRET signal upon cleavage of the FRET probes. FRET has also been used to understand conformational changes in large complexes such as multimeric E3 ligases and theproteasome,frequently using sophisticated single molecule methods. Overall, FRET is a powerful tool to help unravel the intricacies of the complex ubiquitination system.